Mon. Not. R. Astron. Soc. 000, [l]-?? (2013) Printed 6 March 2013 (MN I^TeX style file v2.2) 



Comprehensive abundance analysis of red giants in the 
open clusters NGC 2527, 2682, 2482, 2539, 2335, 2251 and 
2266 

Arumalla B.S. Reddy^*, Sunetra Giridhar^ and David L. Lambert^ 

^Indian Institute of Astrophysics, Bangalore 560034, India 

^W.J. McDonald Observatory, The University of Texas at Austin, Austin, TX 78712 - 0259, USA 



6 March 2013 



ABSTRACT 

We have analyzed high-resolution echelle spectra of red giant members for seven 
open clusters in the Galactic anticentre direction to explore their chemical compo- 
sitions. Cluster membership has been confirmed by radial velocity. The spread in 
temperatures and gravities being very small among the red giants, nearly the same 
stellar lines were employed for all stars thereby reducing the abundance errors: the 
errors of the average abundance for a cluster were generally in the 0.02 to 0.05 dex 
range. Our present sample covers Galactocentric distances of 8.3 to 11.3 kpc and an 
age range of 0.2 to 4.3 Gyr. A careful comparison of our results for the cluster NGC 
2682 (M 67) to other high-resolution abundance studies in the literature shows general 
good agreement for almost all elements in common. 

Key words: - Galaxy: abundances - Galaxy: open clusters and associations ~ stars: 
abundances: general - open clusters: individual: NGC 2251, NGC 2266, NGC 2335, 
NGC 2482, NGC 2527, NGC 2539 and NGC 2682 



1 INTRODUCTION 

This is our second paper reporting abundance measure- 
ments for red giants in open clusters (OCs) lacking de- 
tailed information on their chemical composition. In our 
first paper (Reddy et al. 2012, hereafter paper i) we have 
presented abundance measurements for four open clusters 
(OCs) whose Galactocentric distances (Rgc) lie between 8.3 
and 10.5 kpc, where the Sun's Galactocentric distance is 
taken as 8.0±0.6 kpc (Ghez et al. 2008). Here, we add seven 
OCs six of which have not previously analysed in detail also 
in the Galactic anticentre direction in the disc. Our over- 
all goal is to improve understanding of Galactic chemical 
evolution. 

The abundance variation of chemical elements with 
Galactocentric distance (Rgc) (i.e., the radial metallicity 
gradient) along (and across) the disc and the gradient's tem- 
poral variation over the disc's lifetime put constraints on 
Galactic chemical evolution models which are controlled, in 
large part, by variations of such quantities as the initial mass 
function (IMF), the star formation rate (SFR), gas flows in 
and out of the disc as well as through the disc. Particularly 
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informative diagnostics are the spatial and temporal vari- 
ations of relative abundances which are traditionally refer- 
enced to the iron abundance, i.e., ratios such as O/Fe which 
primarily samples the relative rates of Type II and Type la 
supernovae and Ba/Fe where Ba measures heavy element 
synthesis by AGB stars. 

Abundance variations across the Galaxy, particularly 
measurements of the radial gradient, have been measured 
using a wide variety of objects including Hll regions, hot 
young stars, Cepheid variables, planetary nebulae, red gi- 
ants, and OCs. Various disagreements and inconsistencies 
remain - see, for example, Magrini et al. (2010). OCs either 
through spectroscopy of their red giants or cooler main se- 
quence stars provide not only abundance estimates for many 
elements - essentially, elements sampling all the major pro- 
cesses of stellar nucleosynthesis - but a collection of stars 
with a well determined age, distance and metallicity. More- 
over, it is possible to estimate from an OCs space motion 
and a model of the Galactic gravitational potential the birth- 
place of the OC. 

OCs dissolve over time thanks to tidal forces and feed 
stars into the field. Intracluster abundance variations, espe- 
cially in the a- and the heavy elements provide, in principle, 
a way to chemically tag groups of field stars to their parent 
cluster (Freeman & Bland-Hawthorn, 2002; De Silva et al. 
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2006, 2007). Such studies benefit from tlie largest possible 
sample of OCs with all properties analysed homogeneously 
with high-resolution spectroscopy. 

Addressing abundance variations across the Galactic 
disk and applying chemical tagging requires homogeneous 
and accurate abundance analyses of as large a sample of OCs 
as possible. In this paper, we continue our modest efforts in 
this regard by extending our homogeneous abundance anal- 
ysis for a sample of seven OCs for many elements including 
the s- and r-process elements for which the measurements 
are often lacking in the literature. 

The layout of the present paper is as follows: In Section 
2 we describe the data selection, observations and data re- 
duction, and Section 3 is devoted to the abundance analysis. 
In Section 4 we present our results and compare them with 
the abundances derived from samples of field giants. Finally, 
in Section 5 we give the conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Red giant members of OCs, with the exception of NGC 
2682 (M67), not yet subjected to a comprehensive abun- 
dance analysis using high resolution spectroscopy were se- 
lected from the New catalogue of optically visible open clus- 
ters and candidate^ (Dias et al. 2002). Selection of red gi- 
ants instead of main sequence stars enables the extension of 
our observations to more distant OCs. We have made use of 
the WEBDA[^ database for the selection of suitable red giant 
candidates. 

Target clusters and their properties are shown in the 
Table 1: column 1 represents the cluster name, columns 2 

6 3 the Galactic longitude and latitude in degrees, column 
4 the age, column 5 the photometric estimate of the iron 
abundance, column 6 the Galactocentric distance, column 

7 the distance modulus, column 8 the reddening, column 9 
the reference to the photometric [Fe/H]. All quantities are 
from the database entry except for the photometric [Fe/H] 
abundance and the Galactocentric distance, Rgc, which we 
calculate assuming a distance of the Sun from the Galactic 
centre of 8.0±0.6 kpc (Ghez et al. 2008). 

High-resolution optical spectra of the program stars 
were obtained during the nights of 2011 January 12-14 and 
March 15-18 with the Robert G. TuU echelle coude spec- 
trograph (TuU et al. 1995) on the 2.7-m Harlan J. Smith 
telescope at the McDonald observatory using a 2048 x 2048 
pixel Tektronix charge-coupled device (CCD) as a detector. 
The spectra correspond to a resolving power of > 55,000 (< 
6 km s^^) as measured by the FWHM of Th i lines in com- 
parison spectra. The spectral coverage in a single exposure 
from 4000 A to 5600 A across various orders was complete 
but incomplete from 5600 A to about 9800 A, where the 
inter-order gaps begin to appear. In order to minimize the 
effect of cosmic rays and to acquire a good signal-to-noise 
(S/N) ratio for our spectra, two or three exposures were 
taken each lasting for 20-30 min depending on the magni- 
tude of the star. 

The basic observation procedure and the data reduction 

^ http://www.astro.iag.usp.br/ wilton/ 
^ http://www.univic.ac.at/webda/ 
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Figure 1. Representative spectra of the red giant members of 
the five open clusters described in Table 1 



techniques are same as those listed in our paper i. The echelle 
spectroscopic data is extracted to one-dimensional spectra 
with y-axis as flux and the x-axis as wavelength, using the 
standard spectral reduction software iRAij^ 

Multiple exposure frames from each of the stars were 
median combined to acquire optimal S/N ratio and then 
trimmed to reduce the edge effects during the continuum 
fitting. The final extracted spectra have S/N ratio of about 
100—190 as measured around 6000 A region, while at wave- 
lengths shorter than 5000 A the S/N ratio decreases with 
decreasing wavelength. The continuum fitting was done in- 
teractively by marking continuum regions on an each aper- 
ture which were then fitted by a cubic spline of appropriate 
order for proper continuum fitting. 

We measured the radial velocity (RV) of each star from 
its continuum-fitted spectrum using the cores of weak and 
moderately strong lines of various atomic species. Observed 
radial velocities were converted to heliocentric velocities us- 
ing the rvcorrect routine available in iraf software. Our 
radial velocity measurements are in good agreement with 
the previous radial velocity measurements for the red gi- 
ants in OCs (Mermilliod et al. 2008). But, for NGC 2266, 
no radial velocity measurements from high-resolution spec- 
tra are available in the literature. Recently, Carrera (2012) 
analyzed medium-resolution spectra (R ~ 8000) in the in- 
frared CaT region (~8500 A) for stars in NGC 2266 and 
derived a mean radial velocity of { RV )= — 16±15 km s~^ 

IRAF is a general purpose software system for the reduction 
and analysis of astronomical data distributed by NOAO, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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Table 1. Target clusters and their properties from the hterature. 



Cluster 



(deg.) 



b 

(deg.) 



Age 
(Gyr) 



[Fe/H]phot. 
(dex.) 



Rgc 

(kpc) 



(m-M)v E(B-V) 
(mag.) (mag.) 



[Fe/H], 



NGC 2251 


203.58 


+00.10 


0.27 


-0.20 


9.2 


11.19 


0.19 


Parisi et al. (2005) 


NGC 2266 


187.79 


+10.29 


1.20 


-0.39 


11.3 


12.97 


0.10 


Kaluzny et al. (1991) 


NGC 2335 


223.60 


-01.18 


0.16 


-0.03 


9.1 


11.98 


0.39 


Twarog et al. (1997) 


NGC 2482 


241.62 


+02.03 


0.40 


+0.12 


8.7 


10.93 


0.09 


Twarog et al. (1997) 


NGC 2527 


246.08 


+01.85 


0.45 


-0.09 


8.3 


09.01 


0.04 


Piatti et al. (1995) 


NGC 2539 


233.70 


+11.11 


0.37 


-0.20 


8.9 


10.93 


0.08 


Claria & Lapasset (1986) 


NGC 2682 


215.69 


+31.89 


4.30 


0.00 


8.6 


9.97 


0.06 


Dinescu et al. (1995) 



Table 2. The observed stars 



Cluster 


Star ID 


a(2000.0) 


.5(2000.0) 


V 


B-V 


V-Ks 






S/N 






(hh mm s) 


(O „ ,) 


(mag.) 


(mag.) 


(mag.) 


(mag.) 


(km s^-"-) 


at 6000 A 


NGC 2251 


3 


06:34:51.24 


+08:19:31.90 


10.39 


+1.22 


+2.92 


+0.72 


+26.2+0.2 


150 




33 


06:34:37.07 


+08:21:39.50 


10.39 


+1.21 


+2.91 


+0.70 


+26.3+0.2 


150 


NGC 2266 


73 


06:43:16.69 


+26:57:05.19 


11.06 


+0.99 


+2.46 


+0.69 


-29.7+0.2 


100 


NGC 2335 


11 


07:06:11.42 


-09:56:22.70 


10.89 


+1.13 


+2.76 


+0.67 


-3.21+0.1 


160 


NGC 2482 


9 


07:55:09.09 


-24:22:30.25 


10.27 


+1.11 


+2.38 


+0.61 


+39.00+0.2 


190 


NGC 2527 


10 


08:04:46.97 


-28:07:50.04 


09.49 


+0.95 


+2.17 


+0.59 


+40.7+0.2 


170 




203 


08:05:33.91 


-28:08:58.44 


09.51 


+0.98 


+2.21 


+0.59 


+40.4+0.2 


180 


NGC 2539 


346 


08:10:23.02 


-12:50:43.25 


10.92 


+0.99 


+2.28 


+0.57 


+29.7+0.2 


130 




463 


08:10:42.87 


-12:40:11.80 


10.69 


+1.03 


+2.42 


+0.60 


+28.7+0.2 


135 


NGC 2682 


84 


08:51:12.73 


+ 11:52:42.68 


10.51 


+1.11 


+2.53 


+0.67 


+35.2+0.2 


130 




151 


08:51:26.22 


+ 11:53:52.23 


10.48 


+1.10 


+2.52 


+0.66 


+35.0+0.2 


120 




164 


08:51:29.03 


+ 11:50:33.40 


10.52 


+1.11 


+2.56 


+0.61 


+34.3+0.1 


120 



based on four stars. For this cluster, we selected the star 
with ID 73 whose membership has been confirmed through 
photometry (Kaluzny & Mazur 1991). Our radial velocity 
estimate for this OC would appear to be the first measure- 
ment available in the hterature with negligible dispersion 
whereas Carrera's measurement suffers from large uncer- 
tainty around the mean. Our radial velocity measurement 
for this OC differs from the Carrera's value by —10 km s~^, 
but is well within the quoted uncertainty in Carrera (2012). 

The identification and basic observational data for the 
stars observed in each of the clusters are given in Table 
2 along with the (V-Ks) and (J-Ks) colors from the Two 
Micron All Sky Survey (2MASS) catalogue}^ (Cutri et al. 
2003^ computed radial velocity and S/N of each of the 
spectra extracted at 6000 A for each of the stars. Spectra 
of a representative region are shown in Figure [l] for one star 
from a sample of the five OCs. 



3 ABUNDANCE ANALYSIS 
3.1 Line selection 

Our selection criteria for suitable stellar features from each 
spectra followed the precepts discussed in paper i. The 
equivalent widths (EWs) for each of the selected absorp- 
tion features were measured manually using the routine splot 



contained in iraf by fitting often a Gaussian profile. But 
for lines with significant damping wings a Voigt profile was 
used, for a few lines a direct integration was preferred as a 
best measure of EW. The final linelist includes 250 lines of 
23 elements covering the spectral range 4500 ~ 8850 A. Our 
selection criteria provides, on average across the sample of 
12 stars, a list of 60 Fe i lines with lower excitation poten- 
tials (LEPs) ranging from 0.9 to 5.0 eV, and 12 Fe ii lines 
with LEPs of 2.8 to 3.9 eV and EWs up to 120 mA. The 
line list is essentially that given in Table 3 of Paper I. 

Several lines may be affected by telluric absorption 
lines. The EWs for these lines have been measured only 
when they appear to be Doppler shifted away from the tel- 
luric components, as judged by referring to the Arcturus 
spectrum (Hinkle et al. 2000). 



3.2 Determination of atmospheric parameters 

3.2.1 Photometry 

The initial determination of effective temperature for each 
red giant was derived from dereddenecj^B, V, J and K pho- 
tometry using the empirically calibrated color-temperature 
relations by Alonso et al. (1999). The corresponding errors 
around the (B-V), (V-K) and (J-K) relations are 167 K, 25 



^ http:/ /irsa. ipac.caltech.edu/applications/Gator 
^ Originally published in University of Massachusetts and In- 
frared Processing and Analysis Center (IPAC)/ California Insti- 
tute of Technology. 



^ The adopted interstellar extinctions are (Av, Ax, E(V-K), 
E(J-K))= (3.1, 0.28, 2.75, 0.54)*E(B-V), where E(B-V) is taken 
from WEBDA 
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K and 125 K. Before that the 2MASS magnitudes are 
transformed to standard K magnitude using the relations 
given in Carpenter (2001). The mean difference between the 
two magnitudes is Ks=K+(-0.044±0.003). 

The surface gravities were computed by incorporating 
the known distance to the OCs, effective temperature Teft^ , 
bolometric correction BCv, and the cluster turn-off mass 

into the relation given by (AUende Prieto et al. 1999). 
logg^= log50+log(M*/M0)+4 Iog(r,ff jTes 0) 

+ OA{Vo + BCv) + 2log7v + 0.12 (1) 

with the corresponding luminosity given by 

log{L,/LQ) = -[0.4(Vo + BCv) + 2log tv + 0.12] (2) 

where n is the parallax and Vo is the dereddened Johnson 
V magnitude and the bolometric correction BCv are de- 
rived from estimated effective temperatures and metallicity 
[Fe/H]phot. in Table [l] from the calibration by Alonso et al. 
(1999). We adopt log gQ= 4.44 cm s"^ and Tefi,0= 5777 K. 

A cluster's turn-off mass has been estimated from the 
stellar evolutionary tracks by Yi et al. (2003): turn-off 
masses are 3.3 Mq for NGC 2251, 2.1 M© for NGC 2266, 3.8 
Mq for NGC 2335, 3.0 Mq for NGC 2482, 2.8 Mq for NGC 
2527, 3.1 Mq for NGC 2539 and 1.4 Mq for NGC 2682. Red 
giants are assumed to have the turn-off mass in computing 
the surface gravity. 

Assuming that the various quantities involved in equa- 
tions (1) & (2) are independent of each other and by intro- 
ducing an error of 10% in the stellar mass, an uncertainty of 
3% in Te/ , an uncertainty of 5% in photometric V mag- 
nitude and the bolometric corrections, and an error of 10% 
in the distance (i.e., the parallax), we estimate an error of 
~ 0.11 dex in log with an uncertainty of 0.08 in log L*. 



3.2.2 Spectroscopy 

The spectroscopic abundance analysis was executed using 
the latest version (2010) of the local thermodynamical equi- 
librium (LTE) line analysis and spectrum synthesis code 
MOOG developed by Chris Sneden and originally described 
in Sneden (1973 Model atmospheres were interpolated lin- 
early from the ATLAS9 model atmosphere grid of Castelli & 
Kurucz (2003f] For this purpose, initially we used a model 
with photometrically determined atmospheric parameters. 
These models assume a line-blanketed plane-parallel uni- 
form atmospheres in LTE and hydrostatic equilibrium with 
flux conservation. 

The linelist (Table 13 available online) and basic atomic 
line data (see Table 3 of paper i) give reference solar abun- 
dances that agree well with the published values from As- 
plund et al. (2009). We performed a differential abundance 
analysis relative to the Sun by running the moog in abfind 
mode using the initially estimated stellar parameters (see.. 
Section 3.2.1) and the measured EWs. This abfind driver 



^ The 2MASS survey uses a K-short (Kg) filter whose effective 
wavelength is centred around 2.16 microns in the near-infrared. 
* http://www.as.utexas.edu/ chris/moog.html 
^ http:/ /kurucz. harvard.edu/grids. html 
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Figure 2. The observed trends for the derived abundance from 
the Fe I lines, [Fe I/H], for the star NGC 2682#164 as a function 
of LEP (upper panel), reduced equivalent width (lower panel), 
and the wavelength of each line (lower panel). 



force-fits the individual line abundances to match the com- 
puted EWs to the observed ones, previously measured using 
the splot routine of iraf. 

First, we estimated the microturbulence velocity, ^t, us- 
ing Fe II lines instead of Fe i, since Fe ii lines are less affected 
by the departures from LTE. The is determined by the 
requirement that the abundance from Fe ii line generally 
chosen to have small range in LEP (2.8-3.9 eV) but a good 
range in their EWs, are independent of line's EW or reduced 
EW (log(EW/A). The Te// is estimated by the requirement 
that the Fe abundance from Fe i lines is independent of a 
line's LEP. Finally, the surface gravity, log g, is estimated 
from the constraint that Fe i and Feii lines give the same Fe 
abundance for the derived Tofi and ^t. 

In Figure [2] we show a graphical representation of the 
determination of atmospheric parameters using Fe i line 
abundances as a function of line's LEP, line strength, and the 
wavelength of the lines for a typical star (NGC 2682#164) 
together with a linear fit to each distribution. The vertical 
axis gives the logarithmic iron abundance on the standard 
scale in which log e(H) = 12. The thick dashed blue lines rep- 
resents the mean Fe i abundance, and the thin dashed red 
lines represent linear trends of abundance with the three 
variables. The lack of significant trends of the iron abun- 
dance with LEP (upper panel) and line strength (middle 
panel) supports the validity of our chosen model and the 
adopted ^t. The absence of any trend between the derived 
iron abundance and the wavelength (bottom panel) repre- 
sents a check on our continuum placement. 

A check on the derived spectroscopic atmospheric pa- 
rameters was performed in the same way as described in 
paper i. In short, the uncertainty in the microturbulence 
is estimated from the spread in the derived abundance as 
computed from Fe i/Fe ii, Ti i/Ti ii, V i and Cr i/Cr ii. 
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Table 3. Basic photometric and spectroscopic atmospheric parameters for the stars in each cluster. 



Cluster 


Star ID 




^phot 
cfr V / 




logo 1 . 

& yphot 


rpSpCC 

cfF 


loe Osnrr 


Cspcc 




logfL/Lp)) 






CB-V) 


(V-Kt 


\o ivy 


(cm s ) 




(cm s~^) 


(km sec 


spectroscopy 


photometry 


NGC 2251 


3 


4620 


4642 


4575 


2.21 


4850 


1.90 


1.73 


2.75 


2.36 




33 


4693 


4647 


4624 


2.24 


4850 


2.00 


1.67 


2.65 


2.31 


NGC 2266 


73 


4953 


4910 


4564 


1.80 


4850 


2.60 


1.35 


1.85 


2.67 


NGC 2335 


11 


5372 


5509 


5145 


2.51 


5400 


3.10 


1.68 


1.80 


2.43 


NGC 2482 


9 


4802 


4945 


4756 


2.34 


4850 


2.50 


1.73 


2.10 


2.29 


NGC 2527 


10 


4946 


4994 


4740 


2.80 


5150 


3.10 


1.80 


1.59 


1.84 




203 


4884 


4950 


4734 


2.78 


5050 


2.80 


1.61 


1.85 


1.84 


NGC 2539 


346 


4934 


5002 


4888 


2.63 


5175 


3.10 


1.70 


1.64 


2.05 




463 


4840 


4844 


4794 


2.47 


5050 


2.80 


1.68 


1.90 


2.15 


NGC 2682 


84 


4702 


4683 


4508 


2.36 


4800 


2.60 


1.55 


1.66 


1.85 




151 


4720 


4696 


4552 


2.36 


4700 


2.25 


1.38 


1.98 


1.87 




164 


4702 


4658 


4712 


2.36 


4650 


2.25 


1.48 


1.96 


1.85 



around its minima while the microturbulence is varied over 
the range from to 6 km s~^ (See Figure 2 in our paper i for 
illustration). The uncertainty in the log g value is provided 
by the ionization equilibrium between neutral and ionized 
species of the elements Sc, Ti, V and Cr. Finally, the un- 
certainty in the Tofi is acquired by inspecting the slope of 
the relation between the Fe i abundance and LEPs of the 
lines. Therefore, the typical uncertainties estimated in this 
analysis are 100 K in Te//, 0.25 cm s~^ in log g and 0.20 
km s~^ in ^t. 

The derived stellar parameters for program stars in each 
of the cluster are shown in Table [3] column 1 represents the 
cluster name, column 2 represents the star ID, columns 3 
& 4 represent the photometric Toff using (B-V), (V-K) and 
(J-K) and log g values, columns 5-7 represent the spectro- 
scopic TefT, log g and estimates. Finally, the spectroscopic 
and photometric luminosities ( log(L/L0)) are presented in 
columns 8 & 9. The mean difference in effective tempera- 
tures estimated using (B-V) and (V-K) is —26 ± 64 K and 
using (V-K) and (J-K) is -1-157 ± 121 K. The correspond- 
ing mean differences in ^ , T^f^ and T^ff"^^ with 

respect to spectroscopic T^g'^'^'s are —100 ± 113K, —75 ± 
116K and —232 ± 117K respectively. Mean differences in 
log g and log(L/L0) across the 12 stars are —0.18 ± 0.33 
cm-s~^ and —0.15 ± 0.35 cm-s"'^ respectively. 

The photometric Toff's derived from the calibrations 
based on infrared flux method (Alonso et al. 1999) are sen- 
sitive to the adopted colors, reddening estimates and metal- 
licity. When using the (B-V)— Toff calibration an error of 
0.03 mag on (B-V) and a conservative uncertainty of 20 % in 
reddening translates into a temperature uncertainty of 1.1 
to 1.3 % each. Equivalently, an error of 0.2 dex in [Fe/H]phot 
implies a temperature uncertainty of 1.2 %. When using the 
(V-K)— Toff calibration, we note that an error of 0.03 mag 
on (V-K) and an uncertainty of 20 % in reddening implies 
a temperature uncertainty of 0.1 to 0.7 % each and an error 
of 0.2 dex in [Fe/H]phot implies a temperature uncertainty 
of 0.1 to 0.9 %. Even though the Toff as a function of (J-K) 
has no dependence on metallicity, an error of 0.03 mag on 
(J-K) implies a temperature uncertainty of 2 % while the 
effect of reddening variation is less. As we have shown in 
Paper I, temperatures derived from (V-K) colour might be 
least aflected by photometric uncertainties and we shall use 



these values as the photometric Toff estimates of our giant 
stars. 

With few exceptions, our spectroscopic estimates are 
in good agreement with the photometric ones. The uncer- 
tainty affecting latter is subject to the uncertainties in the 
reddening values. For example, for NGC 2251 the difference 
in T^g and T^^""^ is around 200 K using the reddening 
estimate of E(B-V)= 0.186 as quoted in WEBDA database. 
If we replace this value with the one (E(B-V)= 0.21±0.03) 
given in the recent analysis of Parisi et al. (2005), the differ- 
ence between the two temperatures will be reduced by 50- 
100 K. The corresponding comparison of the spectroscopic 
[Fe/H] with the photometric one in Table 1 also illustrates 
fair agreement: A[Fe/H] = 0.10 (NGC 2251), -0.05 (NGC 
2266), -0.15 (NGC 2335), -0.19 (NGC 2482), -0.01 (NGC 
2527), 0.14 (NGC 2539) and -0.08 (NGC 2682). 

3.3 Abundances and error estimation 

Once the atmospheric parameters were estimated from the 
spectral line measurements as described in previous section, 
the corresponding model atmospheres were used in abfind 
and Synth drives of moog to conduct a complete abundance 
analysis. In most cases, the abundances are derived from the 
measured EWs but a few lines were analysed with synthetic 
spectra. 

Synthetic spectra were computed for species affected by 
hyperfine and isotopic splitting or affected by blends. Llines 
treated by synthesis included: Sc ii (6245 A), Mn i (6013 
A and 6021 A), V i (5727 A), Cu i (5218 A), Rb i (7800 
A), Ba If (5853 A) and Eu ii (6645 A) by including their 
hfs components, and Ce ii (5472 A) and Sm ii (4577 A). 
We adopted the hfs data and isotopic ratios listed in Paper 
I. The hfs data and isotopic ratios for the Eu ii line has 
been kindly provided by C. Sneden whose linelist uses the 
(?/- values provided by Lawler et al. (2001). For the Rb i 
line, the number of hfs components, the accurately known 
wavelengths and the relative strength of each component 
were taken from Lambert & Luck's (1976) analysis of the 
Rb I resonance line in the Arcturus spectrum and a reliable 
log gf value of -f 0.13± 0.04 from Wiese & Martin (1980). A 
fit to this line includes a Si f blend to the blue. As there are 
no reliable measurements of experimental ^/-values for this 
line, we have adopted a solar g f -value of —0.75 to model the 
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Figure 3. Synthetic spectra (continuous lines) and the observed spectrum (dotted line) of NGC 2682 #151 around the Rb I line at 7800 
A. The indicated abundances in the figure are on a logarithmic scale. 



solar Si i line profile from the solar integrated disc spectrum 
(Kurucz et al. 1984). Our solar abundance measurement of 
log £0 = 2.60 dex for the Rb is in close agreement with 
Asplund et al's. (2009) value of log eg = 2.52±0.10 dex. 

A fit to the Rb i line at 7800 A in the observed spectrum 
(dotted line) of NGC 2682 #151 is shown in figure|3]as sohd 
lines with [Rb/H] = 2.35 ± 0.10 dex, based on goodness 
of fit provided by MOOG. A visual inspection of the figure 
shows that the red line ([Rb/H] — 2.35) is a better match 
to the stellar spectrum. 

The largest hfs correction is about —0.11 dex for 
[Sc/Fe], -0.10 dex for [V/Fe], -0.43 dex for [Mn/Fe], -0.39 
dex for [Cu/Fe], -0.12 dex for [Ba/Fe], -0.32 dex for 
[Eu/Fe]. 

Abundance results for the individual stars in each of 
the OCs are presented in Tables [8j-[T2] For each abundance 
based on the analysis of EWs, the abundance and standard 
deviation were calculated from all lines of given species. The 
tables give the [Fe/H] and [X/Fe] for elements considered 
here relative to solar abundances derived from the adopted 
(?/-values. Therefore, errors in the adopted gf values are 
unimportant when providing differential abundances ([X/H] 
or [X/Fe]) provided that the solar and stellar abundances 
depend on the same set of lines. 

Inspection of the Tables [8l- |12| shows that, in general, 
the compositions [X/Fe] of stars in a given cluster are gen- 
erally identical to within the (similar) standard deviations 
computed for an individual star. Exceptions tend to occur 
for species represented by one or a few lines, as expected 
when the uncertainty in measuring equivalent widths is a 
contributor to the total uncertainty. From the spread in the 
abundances for the stars of a given cluster we obtain the 
standard deviation ai in the Tables [S HUI in the column 
headed 'average'. 

We evaluated the sensitivity of the derived abundances 
to the uncertainties in the adopted atmospheric parameters 
by varying each time only one of the parameters by the 
amount corresponding to the typical error. The changes in 
abundances caused by varying atmospheric parameters by 
100 K, 0.25 cm s~^ and 0.2 km s~^ with respect to the 
chosen model atmosphere are summarized in Table |4] We 



Table 4. Sensitivity of abundances to the uncertainties in the model 
parameters for the star with ID 11 in NGC 2335 with Tgff= 5400 K, 



ogg= 3.10 


cm s ^ 


and 5t = 


1.68 km 


s-1. 










TeflltlOO K 


log g± 0.25 


Ct± 0.20 






Species 






<^logg 








Na I 


+0.05/ 


- 0.06 


-0.02/ 


+ 0.02 


-0.04/ 4 


-0.05 


0.04 


Mg I 


+0.03/ 


- 0.03 


-0.01/ 


+ 0.01 


-0.02/ 4 


-0.03 


0.02 


Al I 


+0.03/ 


- 0.04 


-0.01/ 


+ 0.01 


-0.02/ 4 


-0.02 


0.02 


Si I 


+0.01/ 


- 0.01 


+0.03/ 


- 0.02 


-0.02/ 4 


-0.02 


0.02 


Cai 


+0.07/ 


- 0.08 


-0.02/ 


+ 0.02 


-0.09/ 4 


-0.09 


0.07 


Sc I 


+0.09/ 


- 0.09 


0.00/ 


+ 0.01 


0.00/ 4 


0.01 


0.05 


Sc II 


-0.02/ + 0.01 


+0.12/ 


-0.11 


-0.05/ 4 


-0.06 


0.07 


Ti I 


+0.10/ 


- 0.12 


-0.01/ 


0.00 


-0.03/ 4 


-0.03 


0.06 


Ti II 


-0.03/ 


+ 0.01 


+0.11/ 


- 0.12 


-0.07/ 4 


-0.08 


0.08 


V I 


+0.12/ 


- 0.13 


-0.01/ 


0.00 


-0.01/4 


-0.01 


0.05 


Cr I 


+0.08/ 


- 0.10 


-0.01/ 


0.00 


-0.07/ + 0.07 


0.06 


Cr II 


-0.05/ + 0.04 


+0.11/ 


- 0.11 


-0.05/ 4 


-0.06 


0.07 


Mn I 


+0.07/ 


- 0.09 


0.00/ 


- 0.01 


-0.04/ 4 


-0.04 


0.05 


Fe I 


+0.08/ 


- 0.09 


0.00/ 


0.00 


-0.09/ 4 


-0.10 


0.07 


Fe II 


-0.07/ 


+ 0.05 


+0.12/ 


- 0.13 


-0.07/ 4 


-0.07 


0.09 


Co I 


+0.08/ 


- 0.08 


+0.01/ 


- 0.01 


-0.01/ + 0.02 


0.05 


Ni I 


+0.07/ 


- 0.06 


+0.02/ 


- 0.02 


-0.05/ 4 


-0.05 


0.05 


Cu I 


+0.07/ 


- 0.08 


+0.02/ 


- 0.02 


-0.04/ 4 


-0.06 


0.05 


Zn I 


+0.02/ + 0.03 


+0.07/ 


- 0.06 


-0.05/ 4 


-0.07 


0.05 


Rb I 


+0.07/ 


- 0.08 


-0.01/ 


0.00 


-0.01/ 


0.00 


0.04 


Y II 


-0.01/ 


- 0.01 


+0.11/ 


- 0.12 


-0.07/ 4 


-0.09 


0.08 


Zr I 


+0.13/ 


- 0.13 


0.00/ 


- 0.01 


0.00/ 4 


0.01 


0.07 


Zr II 


-0.01/ 


0.00 


+0.10/ 


- 0.11 


-0.01/ + 0.01 


0.06 


Ba II 


+0.01/ 


- 0.03 


+0.08/ 


- 0.10 


-0.17/ + 0.17 


0.11 


La II 


0.00/ 


- 0.02 


+0.10/ 


- 0.12 


-0.02/ 4 


-0.01 


0.06 


Ce II 


-0.01/ 


- 0.01 


+0.10/ 


- 0.11 


-0.01/4 


-0.01 


0.06 


Nd II 


+0.01/ 


- 0.02 


+0.11/ 


- 0.11 


-0.02/ + 0.02 


0.06 


Sm II 


+0.02/ 


- 0.01 


+0.11/ 


- 0.11 


-0.02/ 4 


-0.03 


0.06 


Eu II 


-0.01/ 


0.00 


+0.11/ 


- 0.11 


-0.01/4 


-0.01 


0.06 



quadratically added the three contributors, by taking the 
square root of the sum of the square of individual errors 
associated with uncertainties in temperature, gravity and 
microturbulence, to obtain (T2. The total error atot for each 
of the element is the quadratic sum of cri and a2. The er- 
ror bars in the abundance tables correspond to this total 
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Table 5. Elemental abundance ratios [X/Fe] for elements from Na to Eu for a sample seven OCs from this study 











NOP 9^3q 


NOP 9335 


NnP 9951 




[Na l/Fe] 


+0.32 ± 0.03 


+0.25 + 0.03 


+0.30 + 0.03 


+0.27 + 0.03 


+0.24 + 0.03 


+0.33 + 0.04 


+0.23 + 0.03 


[Mg l/Fe] 


+0.07 ±0.01 


+0.16 + 0.02 


+0.13 + 0.02 


+0.07 + 0.02 


+0.08 + 0.02 


+0.06 + 0.02 


+0.39 + 0.02 


[Al i/Fc] 


+0.05 + 0.02 


+0.09 + 0.01 


+0.07 + 0.02 


0.00 + 0.01 


-0.02 + 0.02 


0.00 + 0.02 


+0.25 ± 0.02 


[Si i/Fe] 


+0.20 + 0.02 


+0.20 + 0.02 


+0.23 + 0.04 


+0.17 + 0.02 


+0.10 + 0.02 


+0.23 + 0.02 


+0.28 + 0.02 


[Ca i/Fe] 


+0.12 + 0.04 


+0.04 + 0.04 


+0.01 + 0.05 


+0.04 + 0.04 


+0.09 + 0.04 


+0.09 + 0.05 


+0.17 + 0.05 


[Sc l/Fe] 


+0.25 ± 0.04 


+0.04 + 0.04 


+0.12 + 0.05 


+0.21 + 0.04 


+0.13 + 0.04 


+0.04 + 0.04 


+0.31 ± 0.05 


[Sc Il/Fe] 


+0.12 ± 0.04 


+0.10 + 0.04 


+0.08 + 0.05 


+0.04 + 0.04 


+ 0.16 


+0.02 


+0.22 ± 0.05 


[Ti i/Fe] 


+0.11 ± 0.04 


-0.01 ± 0.04 


+0.01 ± 0.04 


+0.10 + 0.04 


+0.17 + 0.04 


-0.04 + 0.04 


+0.23 + 0.04 


[Ti ii/Fe] 


+0.08 + 0.04 


+0.01 ± 0.05 


-0.03 + 0.04 


+0.08 + 0.04 


+0.12 + 0.04 


-0.04 ± 0.04 


+0.30 ± 0.04 


[V i/Fe] 


+0.21 ± 0.04 


+0.09 + 0.05 


+0.10 + 0.05 


+0.16 + 0.04 


+0.13 + 0.03 


-0.05 ± 0.05 


+0.20 + 0.06 


[Cr l/Fe] 


+0.10 + 0.03 


+0.05 + 0.03 


+0.10 + 0.04 


+0.08 ± 0.04 


+0.07 + 0.03 


+0.05 ± 0.04 


+0.09 + 0.03 


[Cr Il/Fe] 


+0.06 + 0.04 


+0.08 + 0.05 


+0.09 + 0.05 


+0.08 + 0.04 


0.00 + 0.04 


+0.04 ± 0.05 


+0.02 ± 0.04 


[Mn l/Fe] 


0.00 


-0.08 


-0.11 


+0.01 


-0.01 


-0.13 


-0.01 


[Fe l/H] 


-0.11 ± 0.04 


-0.08 + 0.04 


-0.07 + 0.04 


-0.06 + 0.04 


-0.19+0.04 


-0.10 + 0.05 


-0.45+0.04 


[Fe Il/H] 


-0.09 + 0.05 


-0.08 + 0.05 


-0.07 + 0.05 


-0.07 + 0.05 


-0.17+0.05 


-0.10 + 0.05 


-0.43+0.05 


[Co l/Fe] 


+0.16 + 0.04 


+0.11 + 0.02 


+0.11 + 0.03 


+0.07 + 0.02 


+0.16+0.04 


+0.03 + 0.03 


+0.27 + 0.03 


[Ni i/Fe] 


+0.06 + 0.02 


+0.10 + 0.03 


+0.03 + 0.04 


+0.02 + 0.02 


+0.09+0.03 


+0.04 + 0.03 


+0.09 + 0.03 


[Cu i/Fe] 


-0.14 


0.03 


-0.21 


-0.16 


-0.15 


-0.22 


+0.02 


[Zn l/Fe] 


-0.16 


-0.07 


-0.22 


0.22 


-0.06 


0.13 


0.00 


[Rb l/Fe] 


+0.07 


0.10 


-0.13 


+0.04 


+0.08 


-0.17 


+0.14 


[Y Il/Fe] 


+0.16 + 0.04 


+0.03 + 0.04 


+0.15 + 0.04 


+0.17 + 0.04 


+0.12+0.05 


+0.07 + 0.04 


+0.05 + 0.04 


[Zr l/Fe] 


+0.31 ± 0.05 


-0.07 + 0.05 


+0.11 ± 0.05 


+0.21 ± 0.04 


+0.06+0.04 


0.00 + 0.05 


-0.08+0.05 


[Zr Il/Fe] 


+0.18 + 0.03 


-0.07 + 0.03 


+0.10 + 0.04 


+0.26 + 0.03 


+0.01+0.03 


+0.06 + 0.03 




[Ba Il/Fe] 


+0.08 


-0.16 


+0.09 


+0.10 


+0.25 


+0.11 


-0.13 


[La Il/Fe] 


+0.26 + 0.03 


0.00 ± 0.03 


+0.18 + 0.03 


+0.18 + 0.03 


+0.29+0.03 


+0.02 ± 0.04 


-0.02 ± 0.04 


[Ce Il/Fe] 


+0.24 


-0.02 


+0.11 


+0.20 


+0.29 


+0.07 


-0.07 


[Nd Il/Fe] 


+0.20 + 0.03 


+0.02 + 0.04 


+0.13 + 0.04 


+0.23 + 0.04 


+0.32+0.03 


+0.08 ± 0.04 


+0.13 + 0.04 


[Sm Il/Fe] 


+0.18 


-0.03 


+0.13 + 0.04 


+0.18 


+0.28+0.04 


+0.02 


+0.11 


[Eu Il/Fe] 


+0.10 


+0.08 


+0.07 


+0.19 


+0.07 


+0.04 


+0.40 



Note: Abundances calculated by synthesis are presented in bold numbers. 



error. The final OC mean abundances from this study are 
presented in Table [5] 



4 RESULTS 

The addition of these seven clusters to literature sample 
also supports the widely held impression that there is an 
abundance gradient such that the metallicity [Fe/H] at the 
solar galactocentric distance decreases outwards (Magrini et 
al. 2009, Yong et al. 2012) at about -0.1 dex kpc"^ and at 
the same time our sample has not changed the dispersion of 
metallicities at a given Rgc. 

As in Paper i, we compare our results with those from 
Luck & Heiter (2007) for a large sample of nearby giants 
whose analysis was similar to our, i.e., a differential analysis 
with respect to the Sun. Perhaps, the chief advantage of this 
comparison is that various systematic effects should cancel 
which would enter into a comparison involving our giants 
and field main sequence stars. Two obvious issues are non- 
LTE effects which are, in general, going to be different for 
giants and dwarfs and the use of classical atmospheres which 
simulations of stellar granulation show may in different ways 
be inadequate representation of real stars. 

From Luck & Heiter's Table 4, we calculate mean abun- 
dances [X/Fe] for [Fe/H] from 0.0 to -0.2, the metallicity 
spread of six of our seven clusters. NGC 2266 at [Fe/H] = 
—0.44 is excluded from the sample in part because it may 
belong to the thick rather than the thin disc. One can ar- 



gue based on the membership probabilities (see for example 
Reddy et al. 2006 for definitions and recipes) derived us- 
ing the present day cluster's space motions that this cluster 
has a high probability to belong to the thick disc (Reddy 
& Giridhar, in preparation). Inspection of Table [6] shows 
that the mean [X/Fe] for the six clusters are within 0.10 
dex of Luck & Heiter's results except for Na, Ca, Sc, V, 
Mn, Cu, and Nd and within 0.15 dex for all but Na, Sc, V, 
and Nd. (Elements Zn, Rb, and Sm were not considered by 
Luck & Heiter.) Thus, we conclude that field giants and OCs 
of near-solar metallicity have very similar, if not identical, 
compositions. This was the conclusion reached also in Paper 
I. 

We review our conclusions made in our paper i for Mn i. 
Our [Mn/Fe] ratios for OC giants matches well with Takeda 
(2008) sample of field giants, while it is lower by 0.12 dex 
with respect to the local thin disc field giants of Luck & 
Heiter (2007) . Luck & Heiter's (2007) abundance results for 
field giants are generally confirmed by Takeda et al. (2008) 
for other large sample of field giants. Indeed, Luck & Heiter's 
(2007) [Mn/Fe] abundance ratios for red giants are 0.1 dex 
lower than their sample of field dwarfs (Luck & Heiter 2006). 
Both Luck & Heiter's (2007) and Takeda et al. (2008) have 
employed similar abundance analysis based on the EW mea- 
surement of a large set of lines. Determination of the abun- 
dances of the elements using EW measurements of the lines 
is not always recommended as a good approximation for the 
features affected by hfs components. Moreover, these lines 
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Table 6. Mean elemental abundance ratios, [X/Fe], for Na to 
Eu for the six OCs from this study and the thin disc mean abun- 
dances from Luck & Heiter (2007) in the metallicity of our clus- 
ters (0.0 to —0.2 dex). Abundances calculated by synthesis are 
presented in bold typeface. 



Species 


O C ]Tiea,n 


Thill disc 


[Na l/Fe] 


+0.28 ±0.04 


+0.10 + 0.06 


[Mg l/Fc] 


+0.09 ±0.04 


+0.08 + 0.10 


[Al i/Fe] 


+0.03 + 0.04 


+0.09 + 0.05 


[Si i/Fe] 


+0.19 + 0.05 


+0.12 + 0.04 


[Ca i/Fe] 


+0.07 + 0.04 


-0.04 + 0.05 


[Sc i/Fe] 


+0.13 + 0.09 


-0.08 + 0.06 


[Ti i/Fe] 


+0.06 + 0.08 


0.00 + 0.03 


[V i/Fe] 


+0.11 + 0.09 


-0.09 + 0.07 


[Cr l/Fe] 


+0.07 + 0.02 


+0.01 + 0.05 


[Mn l/Fe] 


-0.05 + 0.06 


+0.06 + 0.07 


[Co l/Fe] 


+0.11 + 0.05 


+0.06 + 0.08 


[Ni i/Fo] 


+0.06 + 0.03 


0.00 + 0.03 


[Cu i/Fo] 


-0.15 + 0.07 


+0.01 + 0.13 


[Rb i/Fe] 


-0.04+0.11 




[Y ii/Fe] 


+0.12 + 0.06 


+0.07 + 0.15 


[Zr i/Fe] 


+0.10 + 0.14 




[Zr Il/Fe] 


+0.09 + 0.12 




[Ba ii/Fe] 


+0.09 + 0.15 


+0.04 + 0.16 


[La ii/Fc] 


+0.15 + 0.12 


+0.05 + 0.09 


[Ce Il/Fe] 


+0.15 + 0.12 


+0.05 + 0.09 


[Nd ii/Fe] 


+0.16 + 0.11 


-0.01 + 0.07 


[Sm Il/Fe] 


+0.13 + 0.11 




[Eu Il/Fe] 


+0.09 + 0.05 


+0.08 + 0.06 



are typically strong in red giant branch stars than their main 
sequence counterparts with similar metallicities. 

Since, the hfs effect desaturates strong absorption lines 
and results in features with larger equivalent widths the 
computed abundances will be overestimated. It will be more 
complicated for the elements having isotopes due to the su- 
perposition of the spectra of the various isotopes of the same 
element. The relative strength and isotopic shift vary from 
element to element. Since, the stellar surface abundance es- 
timates are higly dependent on the line strength, it is im- 
portant to treat strong lines for hfs effects since the hfs cor- 
rection increases with line strength. For weak lines, those 
generally lay on the linear part of the curve of growth, the 
hfs treatment is not so important because they are already 
unsaturated. Hence a spectrum synthesis over the separate 
hfs components is recommended to derive an abundance for 
the lines affected by hfs components. Thus, the number of 
components, the wavelength splittings, isotopic shifts, if any 
present, and the relative strength of each the component 
must be known. 

In our analysis largest hfs corrections are applied to 
[Mn/Fe] (~ -0.43 dex). From the Figure 15 of Takeda (2007) 
for a sample field dwarfs, it is evident that the Mn i lines 
employed for their abundance analysis are all weak (EW < 
50 mA) for which they have noted a negligible hfs correc- 
tions (~ —0.02 dex). Similarly, in their giant sample Takeda 
et al.(2008) have employed a single weak line at 5004.8 A for 
the abundance analysis for which the hfs effects are less im- 
portant. Whereas the Mn i abundances from Luck & Heiter 
(2007), similar to ours, might be based on strong lines for 
which they have not employed hfs treatment, unlike ours. 
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Figure 4. The abundance ratio for [Rb/Zr] as a function of 
[Fe/H] with corresponding error bars. The dots (blue) represent 
the mean cluster abundances derived in this study. The OC NGC 
2266 with a largest [Rb/Zr] ratio (= +0.22) is marked in the 
figure. The empty squares (black) represents the data taken from 
Tomkin & Lambert (1999) for a sample of metal-deficient disc and 
halo stars. The filled square and the filled circle (in red) represents 
the [Rb/Zr] ratios for Arcturus and Sun. The dotted line indicate 
the solar mix of elements. 



This can fairly explain the observed offset in the abundance 
trends of field dwarfs and our OCs red giants with the abun- 
dance trend of field giants from Luck & Heiter (2007). 

The Rb abundance provides an additional diagnostic in- 
formation on the neutron density at the s-process site which 
is controlled by both the neutron source and the mass of the 
parent AGB star responsible for its synthesis. The abun- 
dance of Rb relative to its neighboring s-process elements 
such as Sr, Y and Zr is a useful monitor of neutron den- 
sity at the time of s-processing. Thanks to the branch in 
the s-process path at *^Kr with a /3-decay half life of 10.7 



which at low neutron densities (N„ ^ 10^ cm~^) decays 
to *^Rb, while at higher neutron densities suffers successive 
neutron captures and the flow proceeds along the path from 
**^Kr to ^'^Kr to *^Kr, that ^-decays to stable *''Rb isotope. 
The sensitivity of the Rb elemental abundance to the neu- 
tron density arises then due to the quite different neutron 
capture cross sections for both of its stable isotopes (*^Rb 
and **^Rb) among which the *^Rb being roughly a factor 10 
times more efficient than ^'^Rb. Thus, at low neutron den- 
sities the ratio [Rb/Zr] is about 10 times less than that at 
high neutron densities. Since, the solar Rb abundance is at- 
tributed in almost equal parts to the s-process and r-process, 
the observed stellar Rb abundance directly indicate the neu- 
tron density at the s-process site. Even though the stellar Rb 
I isotopic ratios can not be measured, the Rb i lines used in 
our analysis are weak in all the program stars so the derived 
Rb abundances are not affected by changes in the isotopic 
mixture. We derived the Zr abundances using 5-7 weak lines 
(see the Tables |8- 12 for the number of lines used for each 
star) for which reliable estimates of gf -values are available 
and the usage of the same lines will remove the systematic 
uncertainties. 

A comparison of our observed [Rb/Zr] abundance ratios 
for these OCs (Figure |4| with the recent models of AGB 
nucleosynthesis from Figure 14 in Smith et al. (2000) sug- 
gests that the 3 Mq models provide complete overlap with 
observations except for the [Rb/Zr] ratio in the OC NGC 
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Table 7. Comparison of elemental abundance ratios [X/Fe] for open cluster M 67 using data 
from various sources. 



Species 


This study 


TautvaisieneOO 


Yong05 


PancinolO 


FriellO 


[Na l/Fe] 


+0.25 ±0.02 


+0.20 + 0.00 


+0.30 + 0.05 


+0.08 + 0.09 


+0.13 + 0.06 


[Mg l/Fe] 


+0.16 ±0.02 


+0.09 + 0.00 


+0.16 + 0.05 


+0.27 + 0.04 


+0.05 + 0.06 


[Al i/Fe] 


+0.09 + 0.01 


+0.13 + 0.02 


+0.17 + 0.03 


+0.03 + 0.02 


+0.11 + 0.04 


[Si i/Fe] 


+0.20 + 0.02 


+0.09 + 0.04 


+0.09 + 0.06 


+0.10 + 0.02 


+0.18 + 0.06 


[Ca l/Fe] 


+0.04 + 0.02 


+0.06 + 0.10 


+0.07 + 0.02 


-0.16 + 0.03 


-0.08 + 0.08 


1 _L 1 VJ rii\ 


— u.ui m u.uz 




_l_n 19-1-0 09 
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[Cr i/Fe] 


+0.05 + 0.02 


+0.07 + 0.09 




+0.01 + 0.03 


+0.00 + 0.02 


[Fe i/H] 


-0.08 + 0.02 


0.00 + 0.06 


+0.02 + 0.08 


+0.05 + 0.02 


+0.03 + 0.08 


[Ni l/Fe] 


+0.10 + 0.02 


+0.06 + 0.09 


+0.08 + 0.05 


+0.05 + 0.01 


-0.02 + 0.05 


[Y Il/Fe] 


+0.03 + 0.02 


-0.02 + 0.11 




-0.05 + 0.04 




[Zr l/Fe] 


-0.07 + 0.02 


-0.19 + 0.09 


-0.28 + 0.02 




-0.14 + 0.03 


[Ba Il/Fe] 


-0.16 + 0.06 


+0.08 + 0.00 


-0.02 + 0.05 


+0.25 + 0.02 




[La Il/Fe] 


0.00 + 0.02 


+0.12 + 0.05 


+0.11 + 0.03 


+0.05 + 0.06 




[Ce ii/Fe] 


-0.02 + 0.03 


+0.08 + 0.08 








[Nd ii/Fe] 


+0.02 + 0.03 










[Sm Il/Fe] 


-0.03 + 0.03 


+0.06 + 0.00 








[Eu Il/Fe] 


+0.08 + 0.03 


+0.07 + 0.00 


+0.06 + 0.02 







Note: Abundances calculated by synthesis are presented in bold numbers and the quoted errors 
on them are the sensitivity of their abundance to adopted model atmospheric parameters. An 
internal uncertainty of 0.00 means that the abundance analysis is based on a single line only. 



TautvaisieneOO— Tautvaisiene et al. (2000); Yong05— Yong et al. (2005); PancinolO— Pancino 
et al. (2010); FriellO- Friel et al. (2010); 



2266. For these low mass AGB stars, the dominant neutron 
source is ^^C(oi,nY^O reaction which is responsible for the 
production of s-process elements to observed values. We note 
that NGC 2266 has a higher [Rb/Zr] ratio than the other 
OCs studied here, although our measured Rb abundance is 
based on a single star. Such a high [Rb/Zr] ratio for the OC 
NGC 2266, complimented with sub-solar s-process elemental 
abundances and over-solar r-process and a-elements, reflects 
the addition of Tyepe II SNe to the chemical enrichment of 
its parent gas cloud prior to its formation. 

Our present sample includes the well-studied OC NGC 
2682 (M 67) for which high-resolution spectroscopic abun- 
dance analyses have been reported by Tautvaisiene et al. 
(2000, hereafter TautvaisieneOO), Yong et al. (2005, here- 
after YongOS), Pancino et al. (2010, hereafter PancinolO) 
and Friel et al. (2010, hereafter FriellO). Unfortunately, none 
of the above groups, except TautvaisieneOO sample, have 
stars in common with our sample for a direct star-to-star 
comparison of results, but on the assumption that the in- 
trinsic dispersion in abundances within NGC 2682 is negli- 
gible, we may compare results from the separate analyses. 
For the TautvaisieneOO sample we give in Table [7| the mean 
of the average abundances of three stars in common with our 
study and this mean value does not change if we consider 
the whole sample of stars studied by TautvaisieneOO. Differ- 
ences in [X/Fe] between our results and the four published 
analyses are ±0.15 dex or smaller for almost all elements 
in common: exceptions include Si and Zr for Yong05, Na, 
Ca, and Ba for PancinolO, Ba for TautvaisieneOO, but all 
elements considered by FriellO fall within the ±0.15 limit. 
Further, a direct star-to-star comparison of [X/Fe] abun- 
dance for almost all elements in common between ours and 
TautvaisieneOO fall within ±0.15 dex with the exception of 
[Ba/Fe]. The [Fe/H] estimates for the four published analy- 
ses are consistent but our result is consistently slightly sub- 



solar: [Fe/H] — —0.08 versus 0.00 to +0.05. This comparison 
(Table [7| suggests a quite satisfactory agreement for all el- 
ements given the many different choices which enter into a 
LTE model atmosphere abundance analysis. 

It seems very likely that a concerted effort to establish 
a truly common method of analysis for red giants in OCs 
for use by all the principal players should minimize or even 
eliminate systematic abundance differences. Then, a com- 
bined set of results for OCs across the Galaxy should reveal 
novel insights into the development of abundance differences 
in time and space within the Galactic disc. 



5 CONCLUSIONS 

This study is a part of a project devoted to the measure- 
ment of metallicity and chemical abundances of giant stars 
in OCs covering a broad range of elements, since such mea- 
surements are often lacking in the literature, especially for 
neutron capture elements. We have performed a detailed 
differential abundance analysis of red giants in seven OCs 
located in the Galactic anticentre direction at Galactocen- 
tric distances of Rgc ~ 8.3 to 11.3 kpc. Based on the high 
resolution spectra and standard LTE analysis, we have de- 
rived stellar parameters and abundance ratios for elements 
from Na to Eu. 

The main results of this study are as follows: 

(i) We have derived average [Fe/H] values of — 0.10±0.04 
for NGC 2527, -0.08 ± 0.04 for NGC 2682, -0.07 ± 0.04 
for NGC 2482, -0.06 ±0.04 for NGC 2539, -0.18 ±0.03 for 
NGC 2335, -0.10±0.05 for NGC 2251, and -0.44±0.04 for 
NGC 2266. 

(ii) Our findings confirm the result from Paper i that 
these OC and local disc field giants with [Fe/H] ~ have 
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identical compositions to within the errors of measurements. 
The observed offset of [Mn/Fe] between the Luck & Heiter's 
(2007) sample of field giants and the field dwarfs/OC giants 
might be arising due to the neglect of hfs treatment in their 
abundance analysis of field giants. 
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Table 8. Elemental abundances for stars in the OC NGC 2527. 
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Note: The abundances calculated by synthesis are presented in bold 
numbers. The remaining elemental abundances were calculated using 
line equivalent widths. Numbers in the parentheses indicate the num- 
ber of lines used in calculating the abundance of that element. In this 
analysis we have adopted the hfs data of Prochaska &; McWilliam 
(2000) for Mn I, MucciareUi et al. (2008) for Eu II line, McWilliam 
(1998) for Ba II line, and Allen et al. (2011) for Cu I lines. 
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Table 9. Elemental abundances for stars in the OC NGC 2682. 
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Note:Same as in table m 
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Table 10. Elemental abundances for stars in the OC NGC 2539. 



Species 


star no. 


Q/1 R 
o4D 


star no. 463 


Average 


[Na I/Fe] 


in 97 -1- 


n 

U. 


,01(4) 




97 -1- 


u 


C\A( A\ 


+0.27 + 0.02 


[Mg I/Fe] 


in nfi -t- 


n 

U. 


.05(5) 


-rU. 


\)o in 


u 




+0.07 + 0.02 


[Al I/Fe] 


— u.uu in 


n 

u. 


.02(5) 


-un 

-hU. 




n 


m ( A\ 


0.00 + 0.01 


[Si I/Fe] 


in 1 A -t- 


n 

U. 


.03(12) 


-rU. 


91 -1- 


u 


.UO(^ IZ ) 


+0.17 + 0.02 


[Ca I/Fe] 


n no 4- 


0, 


.04(12) 


1 n 
H-U. 


UD it 


u 


.U4(^10 ) 


+0.04 + 0.03 


[Sc I/Fe] 


n oi -L 
+U.Zi it 


0, 


,03(7) 


1 n 


OI J- 

.21 it 


U 


.U4(^0 ) 


+0.21 ± 0.02 


[Sc II/Fe] 


in -1- 

^U.UO It 


0, 


,06(5) 


-hU. 


Uo It 


u 


.\jOyO ) 


+0.04 ± 0.04 


[Ti I/Fe] 


in 1 n -t- 

^U. lU It 


0, 


,04(15) 


-rU. 


11-1- 
.11 It 


u 


c\A(^ A\ 

.U4(^ 14t ) 


+0.10 + 0.03 


[Ti II/Fe] 


in ns -1- 

^U.UO It 


0, 


,03(7) 


-un 

-hU. 


. uo ^ 


n 


.UOI^U ) 


+0.08 + 0.03 


[V I/Fe] 


in 1 Q -t- 

^U. lO It 


0, 


.05(13) 


-rU. 


1 A -1- 
.14 It 


n 
u 


.U4l^ lO ) 


+0.16 + 0.03 


[Cr I/Fe] 


in 1 n -t- 

^U. lU It 


0, 


,03(10) 


-rU. 


.U / It 


u 


C\A(9.\ 
.U4(^0 ) 


+0.08 + 0.02 


[Cr II/Fe] 


n HQ 4- 

+u.uy It 


0, 


,02(6) 


H-U. 


07 -1- 
.U / it 


u 


HQ 

.Uo(^D ) 


+0.08 + 0.02 


[Mn I/Fe] 


u.uu 






j_n 

-rU 








+0.01 


[Fe I/H] 


n nv -i- 
— U.U/ It 


0, 


,04(51) 


— u. 


.UO it 


u 


.U4(^DD ) 


-0.06 + 0.03 


[Fe II/H] 


n nv -t- 

— U.U/ It 


0, 


,03(13) 


n 

— u. 


07 -1- 
.U / It 


n 
u 


c\A(^ A\ 

.U4l^ 11 j 


-0.07 + 0.02 


[Co I/Fe] 


n 1 o _i_ 
+U.io It 


0, 


,04(4) 


H-U. 


Ol -1- 
.Ul it 


u 


ri/i ( c\\ 

.U4(^D J 


+0.07 + 0.03 


[Ni I/Fe] 


00 it 0. 


,04(17) 


-ru. 


05 i 


Q 


.03(21) 


+0.02 ± 0.02 


[Cu I/Fe] 


-0.18 






-0 


.15 






0.16 


[Zn I/Fe] 


0.23 






-0 


.22 






-0.22 


[Rb I/Fe] 


+0.07 






+0 


.01 






+0.04 


[Y II/Fe] 


+0.15 it 


0, 


.05(6) 


+0. 


.19 it 





.04(5) 


+0.17 + 0.03 


[Zr I/Fe] 


+0.23 it 


0, 


.07(5) 


+0. 


.19 it 





.03(5) 


+0.21 ± 0.04 


[Zr II/Fe] 


+0.28 ± 


0, 


,04(2) 


+0. 


.25 it 


0. 


.03(2) 


+0.26 + 0.02 


[Ba II/Fe] 


+0.13 






+0 


.08 






+0.10 


[La II/Fe] 


+0.18 ± 


0, 


.02(4) 


+0. 


.18 it 





.04(5) 


+0.18 + 0.02 


[Ce II/Fe] 


+0.21 






+0 


.20 






+0.20 


[Nd II/Fe] 


+0.22 + 


0, 


.05(5) 


+0. 


.24 it 


0. 


.03(5) 


+0.23 + 0.03 


[Sm II/Fe] 


+0.16 






+0 


.20 






+0.18 


[Eu II/Fe] 


+0.19 






+0 


.19 






+0.19 



Note:Samc as in table m 
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Table 11. Elemental abundances for stars in the OC NGC 2251. 



Species star no. 3 star no. 33 Average 



[NaT /F^il 




_l_n '^1 -1- C\A( 


in 'i'^ -1- 0"? 
-ru.oo ^ u.uo 


[A/To- T /Ff^l 




_i_n 09-1-0 o'^/'f^^ 


in OR -1- 0"^ 
-pu.uu ^ u.uo 


[AIT 71701 


in n'^ -1- n riA('^\ 


_n 09-1-0 O'^l''^^ 
— u.u^ ^ u.uo^^oy 


n nn -L n 09 
u.uu ^ u.u^ 


[ / J 




_i_n 9f^ -1- c\^(^A\ 
-rU.zo ^ u.uo^^l4:^ 


in 9"^ -1- 09 
-ru.zo ^ u.uz 


rPa T /Ffil 


_i_n 1 A _|_ n f\A(^ fi^ 

-pu. lU ^ U.U^\^ lU J 


_i_n OS -1- o'^/'i 

-ru.uo ^ u.uoi^ioj 


_l_n OQ -1- 09 
-pu.uy ^ u.uz 


[oC 1/ re\ 


1 n no -U n n'i/'cc\ 
-f-U.uy it U.Uol^oj 


— U.Ul it \J.\Joyo) 


in n/1 -U n no 
-1"U.U4 it U.Uz 




_j_Q Q]^ 




+0 02 


[111/ i^ej 


n 01 -1- n oaM 


07 -1- c\A{^ ^\ 
— u.u/ ^ u.un^io^ 


04 -1- 0*^ 

— u.u^ ^ u.uo 


[Ti TT /FpI 

J_ 1 1 1 / 1 c 


u.uu ^ u.u^^^oy 


_n 09-1-0 C\'^(A\ 

— U.U^ ^ U.UtJI^^y 


_n 04 -1- 0*^ 
— U.Ul ^ u.uo 


[V T /Ff>l 
[V 1/ rtij 


n 01 -1- OAM fi^ 

— U.Ul ^ U.U'4^\^1U^ 


n 1 -1- c\'^(^ 

— U.IU ^ U.UOI^lOj 


n 0^1 -1- 09 
— U.UO ^ u.uz 


[Pr T /VcA 


in ns -1- c\^(^ 
-ru.uo ^ u.u•J^^lu^ 


_i_n 09-1-0 o^^/'s^ 


in 0^1 -1- 0*^ 
-[-"U.Utj ^ u.uo 


[Pr TT /FpI 
[V-^l 11/ l^cj 


no -1- o^/'s^ 

u.uu ^ U-UtJl^O^ 


_i_n OS -1- o^/'y^i 


4-0 04 -1- 0*^ 

-pu.u^ ^ u.uo 


[A/Tn T /FpI 

[IVlll 1/ 1^ cj 


— 12 


—0 14 


—0 13 


[Fp T /TTl 


n 1 1 -{- 04/''^fi^ 


n OQ -1- ^^(^7^ 
— u.uy ^ u.uoi^oi J 


n 1 -1- 0*^ 
— U.IU ^ u.uo 


[Fp TT /TTl 

[1^ C 11 / Xlj 


n 1 1 -1- oaM 1 ^ 

— U.ll ^ U.U'4^^^11^ 


n OQ -1- c\A(^ 9^ 

— u.uy ^ u.uii^izj 


n 1 -1- 0*^ 
— u.iu ^ u.uo 


[Pn T /Fol 


_i_n OS -1- of^^'^^^ 


n 01 -1- nAf^\ 

— U.Ul It U.U'4:(^0^ 


in 0*^ -1- 0*^ 
"i-u.uo It u.uo 


[Ni I/Fe] 


+0.05 ±0.04(22) 


+0.04 + 0.04(17) 


+0.04 ±0.03 


[Cu I/Fe] 


-0.19 


-0.27 


-0.22 


[Zn I/Fe] 


-0.11 


-0.16 


0.13 


[Rb I/Fe] 


-0.17 


-0.18 


-0.17 


[Y II/Fe] 


+0.05 ±0.03(5) 


+0.10 + 0.05(5) 


+0.07 ±0.03 


[Zr I/Fe] 


+0.03 + 0.05(5) 


-0.02 + 0.03(4) 


0.00 ±0.03 


[Zr II/Fe] 


+0.06 + 0.01(2) 


+0.06 + 0.02(2) 


+0.06 + 0.01 


[Ba II/Fe] 


+0.12 


+0.10 


+0.11 


[La II/Fe] 


+0.05 + 0.04(5) 


0.00 + 0.03(3) 


+0.02 ± 0.02 


[Ce II/Fe] 


+0.07 


+0.08 


+0.07 


[Nd II/Fe] 


+0.08 + 0.03(9) 


+0.09 + 0.04(8) 


+0.08 + 0.02 


[Sm II/Fe] 


+0.03 


+0.01 


+0.02 


[Eu II/Fe] 


+0.05 


+0.04 


+0.04 



Note:Same as in table m 
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NGC 2527, NGC 2682, NGC 2482, NGC 2539, NGC 2335, NGC 2251, NGC 2266 



Table 12. Elemental abundances for stars in the OCs NGC 2482, 
NGC 2335 and NGC 2266. 



Species star no. 2482#9 star no. 2335#11 star no. 2266#73 



riVa T /■pc.i 

1 i\ d 1 / Jr c 1 


-un 

-j-U. 






-un 9A -1- n n9/'fi^ 


in 9"^ -1- n r\'^( ^\ 


[Wig i/i^ej 


-rU. 


.10 in u. 




_i_n ns -1- n ^9/'fi^ 


in "^Q -1- n C\0[^\ 
-j-u.oy m u.uzi^oj 


FAIT /TTol 


-rU. 
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.u 1 m vj. 
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[■^i T /Pd 
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-un 

-j-U. 




.U'i^ lO ^ 


_i_n 1 n -1- n n9/'i 


_i_n OQ _i_ n n9('i 


fPo T /Fol 
yK^di, 1/ rtij 


-rU. 


.ui m u. 


.UOI^ 10 ) 


_i_n HQ -1- n fizi/'i 


-ru.ii mu.uoi^iDj 


[oc 1/rej 


-rU. 


19-1-0 


.UOl^D ) 


_i_n 1 -1- n nAf'^\ 




rc^ TT /Fol 


-rU. 








in 99 -1- n 
-j-u.zz m u.uoi^oj 


1 _L i 1 / P 


4-0 
-j-U. 


ni -1- n 

.Ul u. 


HAM 


_i_n 1 7 -1- n c\A(^ a\ 


_i_n 9"^ -1- n nA( i fi^i 


1 _L i 11 /PCI 


n 

u. 


.uo ^ u. 




4_ n 1 9 -1- n fizi/'i 


in "^n -1- n nA( 7\ 


rv T / Pol 

[V 1/ rcj 


-j-U. 


1 n -1- 

. lU ^ u. 


.uo^^ lu ^ 


_i n 1 -1- n c\A(^ c\\ 

"T~U.lO ^ U-UH^lUj 


in 90-1-0 Ofi/' 1 7'! 


rPv T /TTril 

[^r i/rtij 




1 n -1- n 

.lu m u. 


.U^^:!^ IZ ^ 


1 n 07 -1- n ^'^/'1^^ 


in no -1- n n'^/'i9~i 


rpv TT 

11/ rej 










in 09 -1- fl ^zl/'Q^ 


[iviii 1/ 1* ej 


n 






n ni 


n m 
— u.ux 


fPc. T /TTl 

[pe 1/ nj 


— u. 






n 1 Q-i-n nA(^Q\ 


— u.'yiOiiiu.u4t(^oo j 


[Fp TT / TTl 

IP c 11 / HI 


n 






n 1 7-1-0 f)'i( ^ 9^ 

— U. 1 / ^U.UO^^ IZ J 


n 4*^-1-0 O'^/'l 9^1 
— U.^lO^U.UOI^-LZ } 


[Co 1/Fe] 


+0, 


,11 ± 0, 


,03(7) 


+0.16±0.03(3) 


+0.27 ±0.03(5) 


[Ni I/Fe] 


+0, 


,03 ±0, 


.04(17) 


+0.09±0.03(21) 


+0.09 ±0.03(20) 


[Cu I/Fe] 





.21 




-0.15 


+0.02 


[Zn I/Fe] 


-0 


.22 




-0.06 


0.00 


[Rb I/Fe] 


-0 


.13 




+0.08 


+0.14 


[Y II/Fe] 


+0, 


,15 ±0, 


,04(5) 


+0.12±0.05(5) 


+0.05 ± 0.04(4) 


[Zr I/Fe] 


+0, 


,11 ± 0, 


,05(5) 


+0.06±0.04(2) 


-0.08±0.05(5) 


[Zr II/Fe] 


+0, 


,10 ±0, 


,04(3) 


+0.01 ±0.02(2) 




[Ba II/Fe] 


+0 


.09 




+0.25 


0.13 


[La II/Fe] 


+0, 


,18 ±0, 


,03(5) 


+0.29+0.03(4) 


-0.02 ± 0.04(6) 


[Ce II/Fe] 


+0 


.11 




+0.29 


-0.07 


[Nd II/Fe] 


+0, 


,13 ±0, 


,04(7) 


+0.32+0.03(8) 


+0.13 + 0.04(4) 


[Sm II/Fe] 


+0, 


,13 ±0, 


.04(4) 


+0.28+0.03(3) 


+0.11 


[Eu II/Fe] 


+0 


.07 




+0.07 


+0.40 



Note:Same as in table H 
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